Search for: All records

Abstract The LMC’s stellar bar is offset from the outer disk center, tilted from the disk plane, and does not drive gas inflows. These properties are atypical of bars in gas-rich galaxies, yet the LMC bar’s strength and radius are similar to typical barred galaxies. UsingN-body hydrodynamic simulations, we show that the LMC’s unusual bar is explainable if there was a recent collision (impact parameter ≈2 kpc) between the LMC and SMC. Pre-collision, the simulated bar is centered and coplanar. Post-collision, the simulated bar is offset (≈1.5 kpc) and tilted (≈8 $\stackrel{\text{°}}{.}$6). The simulated bar offset reduces with time, and comparing with the observed offset (≈0.8 kpc) suggests the timing of the true collision to be 150–200 Myr ago. Then, 150 Myr post-collision, the LMC’s bar is centered with its dark matter (DM) halo, whereas the outer disk center is separated from the DM center by ≈1 kpc. The SMC collision produces a tilted-ring structure for the simulated LMC, consistent with observations. Post-collision, the simulated LMC bar’s pattern speed decreases by a factor of 2. We also provide a generalizable framework to quantitatively compare the LMC’s central gas distribution in different LMC–SMC interaction scenarios. We demonstrate that the SMC’s torques on the LMC’s bar during the collision are sufficient to explain the observed bar tilt, provided the SMC’s total mass within 2 kpc was (0.8–2.4) × 10⁹M_⊙. Therefore, the LMC bar’s tilt constrains the SMC’s pre-collision DM profile, and requires the SMC to be a DM-dominated galaxy. 

Abstract The first infall of the LMC into the Milky Way (MW) represents a large and recent disruption to the MW circumgalactic medium (CGM). In this work, we use idealized, hydrodynamical simulations of an MW-like CGM embedded in a dark matter halo with an infalling LMC-like satellite initialized with its own CGM to understand how the encounter is shaping the global physical and kinematic properties of the MW CGM. First, we find that the LMC drives order-unity enhancements in MW CGM density, temperature, and pressure due to a $\mathcal{M}\approx 2$shock from the supersonic CGM–CGM collision. The resulting shock front extends from the LMC to beyond ∼R_200,MW, amplifying column densities, X-ray brightness, thermal Sunyaev–Zeldovich distortion, and potentially synchrotron emission from cosmic rays over large angular scales across the southern hemisphere. Second, the MW’s reflex motion relative to its outer halo induces a dipole in CGM radial velocities, withv_R ± 30–50 km s⁻¹atR > 50 kpc in the northern and southern hemispheres, respectively, consistent with measurements in the stellar halo. Finally, ram pressure strips most of the LMC’s CGM, leaving ∼10⁸⁻⁹M_⊙warm ionized gas along the past orbit of the LMC, moving at high radial and/or tangential velocities ∼50–100 kpc from the MW. Massive satellites like the LMC leave their mark on the CGM structure of their host galaxies, and signatures of such interactions may be observable in key all-sky tracers of the MW CGM and those of other massive galaxies. 

Abstract The large-scale morphology of Milky Way (MW)–mass dark matter halos is shaped by two key processes: filamentary accretion from the cosmic web and interactions with massive satellites. Disentangling their contributions is essential for understanding galaxy evolution and constructing accurate mass models of the MW. We analyze the time-dependent structure of MW-mass halos from zoomed cosmological-hydrodynamical simulations by decomposing their mass distribution into spherical harmonic expansions. We find that the dipole and quadrupole moments dominate the gravitational power spectrum, encoding key information about the halo’s shape and its interaction with the cosmic environment. While the dipole reflects transient perturbations from infalling satellites and damps on dynamical timescales, the quadrupole—linked to the halo’s triaxiality—is a persistent feature. We show that the quadrupole’s orientation aligns with the largest filaments, imprinting a long-lived memory on the halo’s morphology even in its inner regions (∼30 kpc). At the virial radius, the quadrupole distortion can reach 1–2 times the spherical density, highlighting the importance of environment in shaping MW-mass halos. Using multichannel singular spectrum analysis, we successfully disentangle the effects of satellite mergers and filamentary accretion on quadrupole. We find that, compared to isolated MW–LMC simulations that typically use a spherical halo, the LMC-mass satellite induces a quadrupolar response that is an order of magnitude larger in our cosmological halo. This highlights the need for models that incorporate the MW’s asymmetry and time evolution, with direct consequences for observable structures such as disk warps, the LMC-induced wake, and stellar tracers—particularly in the era of precision astrometry. 

Abstract The infall of the Large Magellanic Cloud (LMC) into the Milky Way’s halo impacts the distribution of stars and dark matter (DM) in our Galaxy. Mapping the observational consequences of this encounter can inform us about the properties of both galaxies, details of their interaction, and possibly distinguish between different DM models.N-body simulations predict a localized overdensity trailing the LMC’s orbit both in baryonic and DM, known as the wake. We collected wide-field, deep near-infrared, and optical photometry using VIRCAM and DECam across four fields along the expected wake, covering the sky region expected to span most of its predicted density contrast. We identify over 400 stars comprising two different tracers, near main-sequence turnoff stars and red giants, which map the halo between 60 and 100 kpc, deriving stellar halo densities as a function of sky position and Galactocentric radius. We detect (1) a break in the halo radial density profile at 70 kpc not seen in northern halo studies and (2) a clear halo overdensity starting also at 70 kpc, with density contrast increasing steadily toward the expected current location of the wake. If this overdensity is the LMC wake, its peak density contrast is as pronounced as the most massive LMC model considered. Contamination from unidentified substructures may bias our wake detections, so wider-area surveys with similar depth are needed for confirmation. 

Abstract We present new measurements of the 2D geometry of the LMC’s stellar bar with precise astrometric observations of red clump stars in Gaia DR3. We develop a novel solution to tackle crowding-induced incompleteness in Gaia data sets with the Gaia BP-RP color excess. Utilizing the color excess information, we derive a 2D completeness map of the LMC’s disk. We find that incompleteness biases the bar measurements and induces large uncertainties. With the completeness-corrected 2D red clump map, we precisely measure the LMC bar’s properties using Fourier decomposition. The bar radius (semimajor axis) is $R_{\mathrm{bar}}={2.13}_{-0.04}^{+0.03}$kpc, and its position angle is 121.°26 ± 0.°21. The bar’s strength as quantified by the Fourier bi-symmetric amplitude isS_bar= 0.27, indicating that the LMC has a significant bar perturbation. We find the bar has an axis ratio of 0.54 ± 0.03, and is offset with respect to the center of the outer disk isophote atR≈ 5 kpc by 0.76 ± 0.01 kpc. These LMC bar properties agree with a hydrodynamic model where the SMC has undergone a recent direct collision with the LMC. We compare the LMC’s bar properties with other barred galaxies in the local Universe, and discover that the LMC is similar to other barred galaxies in terms of bar-galaxy scaling relations. We discuss how our completeness correction framework can be applied to other systems in the Local Group. 

Abstract The merger timescales of isolated low-mass pairs (10⁸<M_*< 5 × 10⁹M_⊙) on cosmologically motivated orbits have not yet been studied in detail, though isolated high-mass pairs (5 × 10⁹<M_*< 10¹¹M_⊙) have been studied extensively. It is common to apply the same separation criteria and expected merger timescales of high-mass pairs to low-mass systems, however, it is unclear if their merger timescales are similar, or if they evolve similarly with redshift. We use the Illustris TNG100 simulation to quantify the merger timescales of isolated low-mass and high-mass major pairs as a function of cosmic time, and explore how different selection criteria impact the mass and redshift dependence of merger timescales. In particular, we present a physically motivated framework for selecting pairs via a scaled separation criterion, wherein pair separations are scaled by the virial radius of the primary’s Friends-of-Friends (FoF) group halo (r_sep< 1R_vir). Applying these scaled separation criteria yields equivalent merger timescales for both mass scales at all redshifts. Alternatively, static physical separation selections applied equivalently to all galaxy pairs at all redshifts lead to a difference in merger rate of up to ∼1 Gyr between low- and high-mass pairs, particularly forr_sep< 150 kpc. As a result, applying the same merger timescales to physical-separation-selected pairs will lead to a bias that systematically overpredicts low-mass galaxy merger rates. 

Abstract Stellar streams in the Milky Way are promising detectors of low-mass dark matter (DM) subhalos predicted by ΛCDM. Passing subhalos induce perturbations in streams that indicate the presence of the subhalos. Understanding how known DM-dominated satellites impact streams is a crucial step toward using stream perturbations to constrain the properties of dark perturbers. Here, we cross-match a Gaia Early Data Release 3 and SEGUE member catalog of the Cetus-Palca stream (CPS) with H3 for additional radial velocity measurements and fit the orbit of the CPS using this six-dimensional (6D) data. We demonstrate for the first time that the ultra-faint dwarf Segue 2 had a recent (77 ± 5 Myr ago) close flyby (within the stream's 2σwidth) with the CPS. This interaction enables constraints on Segue 2’s mass and density profile at larger radii ( $O\left(1\right)$kpc) than are probed by its stars ( $O\left(10\right)$pc). While Segue 2 is not expected to strongly affect the portion of the stream covered by our 6D data, we predict that if Segue 2’s mass within ∼ 6 kpc is 5 × 10⁹M_⊙, the CPS's velocity dispersion will be ∼ 40 km s⁻¹larger atϕ₁ > 20° than atϕ₁ < 0°. If no such heating is detected, Segue 2’s mass cannot exceed 10⁹M_⊙within ∼ 6 kpc. The proper motion distribution of the CPS near the impact site is mildly sensitive to the shape of Segue 2’s density profile. This study presents a critical test for frameworks designed to constrain properties of dark subhalos from stream perturbations. 

ABSTRACT We obtain a quantitative star formation history (SFH) of a shell-like structure (‘shell’) located in the northeastern part of the Small Magellanic Cloud (SMC). We use the Survey of the MAgellanic Stellar History to derive colour–magnitude diagrams (CMDs), reaching below the oldest main-sequence turnoff, from which we compute the SFHs with CMD-fitting techniques. We present, for the first time, a novel technique that uses red clump (RC) stars from the CMDs to assess and account for the SMC’s line-of-sight depth effect present during the SFH derivation. We find that accounting for this effect recovers a more accurate SFH. We quantify an $$\sim$$7 kpc line-of-sight depth present in the CMDs, in good agreement with depth estimates from RC stars in the northeastern SMC. By isolating the stellar content of the northeastern shell and incorporating the line-of-sight depth into our calculations, we obtain an unprecedentedly detailed SFH. We find that the northeastern shell is primarily composed of stars younger than $$\sim$$500 Myr, with significant star formation enhancements around $$\sim$$250 and $$\sim$$450 Myr. These young stars are the main contributors to the shell’s structure. We show synchronicity between the northeastern shell’s SFH with the Large Magellanic Cloud’s (LMC) northern arm, which we attribute to the interaction history of the SMC with the LMC and the Milky Way (MW) over the past $$\sim$$500 Myr. Our results highlight the complex interplay of ram pressure stripping and the influence of the MW’s circumgalactic medium in shaping the SMC’s northeastern shell. 

Abstract We analyze four epochs of Hubble Space Telescope imaging over 18 yr for the Draco dwarf spheroidal galaxy. We measure precise proper motions for hundreds of stars and combine these with existing line-of-sight (LOS) velocities. This provides the first radially resolved 3D velocity dispersion profiles for any dwarf galaxy. These constrain the intrinsic velocity anisotropy and resolve the mass–anisotropy degeneracy. We solve the Jeans equations in oblate axisymmetric geometry to infer the mass profile. We find the velocity dispersion to be radially anisotropic along the symmetry axis and tangentially anisotropic in the equatorial plane, with a globally averaged value $\overline{{\beta }_{\mathrm{B}}}=-{0.20}_{-0.53}^{+0.28}$, (where 1 – ${\beta }_{\mathrm{B}}\equiv 〈v_{\tan }^2〉/〈v_{\mathrm{rad}}^2〉$in 3D). The logarithmic dark matter (DM) density slope over the observed radial range, Γ_dark, is $-{0.83}_{-0.37}^{+0.32}$, consistent with the inner cusp predicted in ΛCDM cosmology. As expected given Draco’s low mass and ancient star formation history, it does not appear to have been dissolved by baryonic processes. We rule out cores larger than 487, 717, and 942 pc at 1σ, 2σ, and 3σconfidence, respectively, thus imposing important constraints on the self-interacting DM cross section. Spherical models yield biased estimates for both the velocity anisotropy and the inferred slope. The circular velocity at our outermost data point (900 pc) is ${24.19}_{-2.97}^{+6.31}\mathrm{km}{\mathrm{s}}^{-1}$. We infer a dynamical distance of ${75.37}_{-4.00}^{+4.73}$kpc and show that Draco has a modest LOS rotation, with $\left(v/\sigma \right)=0.22\pm 0.09$. Our results provide a new stringent test of the so-called “cusp–core” problem that can be readily extended to other dwarfs. 

Abstract The interaction between the supersonic motion of the Large Magellanic Cloud (LMC) and the circumgalactic medium (CGM) is expected to result in a bow shock that leads the LMC’s gaseous disk. In this letter, we use hydrodynamic simulations of the LMC’s recent infall to predict the extent of this shock and its effect on the Milky Way’s (MW) CGM. The simulations clearly predict the existence of an asymmetric shock with a present-day standoff radius of ∼6.7 kpc and a transverse diameter of ∼30 kpc. Over the past 500 Myr, ∼8% of the MW’s CGM in the southern hemisphere should have interacted with the shock front. This interaction may have had the effect of smoothing over inhomogeneities and increasing mixing in the MW CGM. We find observational evidence of the existence of the bow shock in recent Hαmaps of the LMC, providing a potential explanation for the envelope of ionized gas surrounding the LMC. Furthermore, the interaction of the bow shock with the MW CGM may also explain the observations of ionized gas surrounding the Magellanic Stream. Using recent orbital histories of MW satellites, we find that many satellites have likely interacted with the LMC shock. Additionally, the dwarf galaxy Ret2 is currently sitting inside the shock, which may impact the interpretation of the reported gamma-ray excess in Ret2. This work highlights how bow shocks associated with infalling satellites are an underexplored yet potentially very important dynamical mixing process in the circumgalactic and intracluster media.